Nuclear spin pumping and electron spin susceptibilities 
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In this work we present a new formalism to evaluate the nuclear spin dynamics driven by hyperfine 
interaction with non-equilibrium electron spins. To describe the dynamics up to second order in the 
hyperfine coupling, it suffices to evaluate the susceptibility and fluctuations of the electron spin. Our 
approach does not rely on a separation of electronic energy scales or the speciflc choice of electronic 
basis states, thereby overcoming practical problems which may arise in certain limits when using a 
more traditional formalism based on rate equations. 



INTRODUCTION 

In recent years, considerable theoretical and experi- 
mental work is aimed at the controlled manipulation of 
electron spins in nanoscale solid state devices This 
research is motivated by actual applications, such as in 
digital information storage and read-out [2, Q , but also 
by the possibility of using the spin of electrons as com- 
putational units (qubits) in a quantum computer [ij . 

One of the mechanisms influencing the electron spin 
dynamics in these nano-devices is the hyperfine inter- 
action between the electron spin and the nuclear spins 
of the device's constituent material. For spin qubits 
hosted in semiconductor quantum dots Q, hyperfine in- 
teraction has been identified as the main source of deco- 
herence, causing the spin coherence time to be in the 
ns range 0-11] ■ Much recent experimental and theo- 
retical work is aimed at suppressing this hyperfine in- 
duced decoherence 0, • In other semiconductor nano- 
structures [ill, El] — possibly also metallic structures [l3] 
— hyperfine interaction can even dominate the electronic 
transport properties and spin dynamics. Understand- 
ing the role of hyperfine interaction in spintronic devices 
therefore is crucial in the development of the field. 

Traditionally, the interplay between electronic and nu- 
clear spins is treated in the convenient framework of rate 
equations , [Tij . The rates of hyperfine transitions flip- 
ping nuclear spins 'up' and 'down' are separately cal- 
culated using Fermi's golden rule, and the balancing of 
these rate yields the net nuclear spin pumping in the 
system. Although this approach works very well in many 
cases , , it can become cumbersome when the elec- 
tron spin dynamics are more complicated. Strong dy- 
namical nuclear spin effects have been observed in GaAs 
double quantum dots in the spin blockade regime un- 
der conditions of electron spin resonance [l3|- Nuclear 
spin flips in this setup are due to second order processes, 
and some of the transitions have a vanishingly small en- 
ergy difference between initial and virtual state, which 
cannot be dealt with in a standard Fermi golden rule 
approach [l^]- The situation is even worse for systems 
with strong spin-orbit coupling. In InAs nanowire double 



quantum dots in the spin blockade regime, signatures of 
strong dynamic polarization have been observed too (20| . 
In this case, the presence of too many comparable energy 
scales makes it impossible to choose a unique set of ba- 
sis states to describe the electron dynamics in 2l|. It 
is clear that a basis-independent description of the cou- 
pled electron-nuclear spin dynamics, not involving any 
separate transition rates, is highly desirable. 

In this work we present an alternative method to eval- 
uate the spin dynamics of electrons and nuclei coupled 
to each other by hyperflne interaction. We show that, in 
order to describe the nuclear spin dynamics up to second 
order in the hyperfine coupling, it suffices to evaluate the 
ffuctuations and susceptibility of the electron spins in the 
system, which can be done using linear response theory. 
Our approach does not rely on the calculation of sepa- 
rate transition rates with Fermi's golden rule, nor on the 
choice of electronic basis states to work in. As a result, 
our formalism may be applicable in cases where Fermi's 
golden rule cannot be used, and thus provides a useful 
alternative. 

Since a Fermi's golden rule (FGR) calculation often 
has to be adapted to specific circumstances 3-121, it is 
difficult to pick a single implementation of it and make 
a good comparison to our approach. Therefore we will 
illustrate our ideas with two example systems: (i) We will 
show how in a simplest toy system our formalism and a 
smart implementation of FGR produce identical results, 
(ii) We will investigate a more complicated system, which 
cannot be dealt with in an FGR approach, and show how 
our formalism straightforwardly produces an equation for 
the hyperfine driven nuclear spin dynamics. 



MAIN RESULT 

Let us start with presenting the main result of our 
work: an equation for the hyperfine driven dynamics of 
the expectation value of a nuclear spin coupled to (many) 
electronic spins. Under conditions which we will specify 
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below, this equation takes the closed form 



dt 



^abc oh 
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_pa 



(1) 



where {a, 5, c} G {x,y,z}, and we use the convention 
that over repeated indices still has to be summed. The 
operator K is the nuclear spin operator for the nucleus 
under consideration, and e"*"^ is the permutation tensor. 
The vectors S and P, and the matrix Q are defined as 



5'^ = Afo(5"(r„,t)) 



(2) 



A being the material-specific hyperfine coupling energy, 
1 /wo the density of nuclei, K the total nuclear spin, r„ the 
position of the nucleus, and S the electron spin density 
operator, S(r,i) = \'ijj]^{r,t)<Tai3^p{r,t). Note that we 
have set for transparency ?i = 1, or, in other words, we 
express all energies in terms of corresponding frequencies. 
We further assumed for simplicity only one type of spin- 
carrying nuclei present, having a constant density. The 
symbols x R represent correlation functions of the 
local electronic spin density. 



([^"(r„,i),^^(r„,i')] + )cd<' 

(3) 



x''^ = -i / ([^'^(r„,i),5^(r„,i')]->crfi' 



where the square brackets denote the (anti)commutator 
of two operators, i.e., [A, B]± = AB ± BA, and the sub- 
script c means that we have to use the 'connected' part 
of the expectation value, i.e., remove the contribution of 
the averages {AB)^ = {AB) - {A){B). 

Note that, since hyperfine interaction works two-way, 
the electronic dynamics, and thereby S, Q and P, can in 
turn also be affected by the state of the nuclear spins in 
the system. This creates a feedback mechanism in Eq. 
p^, which makes the equation non-linear. However, in all 
practical implementations, one can make use of the large 
difference in time scales of the electronic and nuclear spin 
dynamics. This allows to treat the ensemble of nuclear 
spins effectively as a static 'classical' magnetic field when 
evaluating the electronic correlators. 

The first term in Eq. ([1]) is first order in the hyperfine 
coupling A and gives rise to a precession of the nuclear 
spin around the direction of the time-averaged local elec- 
tron spin polarization. The precession frequency depends 
on the magnitude of this polarization |tio(S(r„, as 
well as on the strength A of the hyperfine coupling. 

The other two terms are second order in A, and de- 
scribe the effect on the nucleus of fluctuations of electron 
spin. The correlation functions i?"'' in the matrix Q are 
ordered as the 'classical' noise of electron spin S(r„, t). It 



persists even if the quantum operator of spin is replaced 
with a classical fluctuating field affecting the nuclear spin. 
Since such fluctuations cause random rotation of this spin 
and therefore isotropization of its density matrix, the ef- 
fect of is relaxation of the nuclear spin. The cor- 
relators in the vector P have the same structure as 
a Kubo formula and give the stationary susceptibility, 
that is, the response of Sa in the point r„ to a con- 
stant magnetic field in the 6-direction, concentrated in 
the same point. One can interpret the x"'^ ^is the 'quan- 
tum' noise of electron spin: the part of the fluctuations in 
S(r„. t) which solely exists due to the non-commutativity 
of the spin operators. It is worth to note that pumping 
of the nuclear spin is produced by the asymmetric part 
of the tensor x""^ only. The Onsager relations state that 
susceptibilities are symmetric for an equilibrium system, 
X"^ = X^°- Therefore, pumping can only take place if the 
electron part of the system is driven out of equilibrium 
by some external agent. 

In order to find an explicit expression for (i(K) /dt as 
given by Eq. ([T]), one needs to find the elements of S, Q 
and P. The first step is to find the steady-state electron 
spin density (S(r„)), which immediately yields the vec- 
tor S. The elements of Q and P correspond to the local 
fluctuations and susceptibility of this electron spin den- 
sity, and can be evaluated using linear response theory, 
as we will illustrate below. Let us emphasize here that 
none of the steps in such a calculation involves the eval- 
uation of any separate transition rate between different 
electronic-nuclear levels or depends on the specific choice 
of electronic basis states. 



DERIVATION 

We will now show how to derive Eq. ([T]) from a second 
order perturbation theory in the hyperfine interaction. 
We start from an effective equation of motion for the re- 
duced density matrix jjf' of the nuclear spins which reads 
up to second order in the perturbation Hhf(t) (see the 
Appendix for a derivation) 



dp" ^ 



[HMit)dHh{it'),p'o®p''U-dt'Y 



(4) 



where is the electronic part of the (quasi-)stationary 
density matrix described by the unperturbed Hamilto- 
nian Hq- This Hamiltonian includes all processes and 
degrees of freedom which are relevant for the electron 
spin dynamics, such as (e.g. in a transport setup) all leads 
and tunnel couplings to these leads. The trace Tre{. . . } 
denotes a trace over the electronic degrees of freedom, 
where the horizontal bar indicates that from the result- 
ing correlation functions only the fully connected contri- 
bution should be taken into account (see Appendix). 
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The perturbation Hh{{t) accounts for the hyperfine in- 
teraction which couples the nuclear spins to the local 
electron spin density. Assuming that the electrons of in- 
terest have s-type orbitals, we use 

HM{t)=AvoY,S''irn,t)k^{t), (5) 

n 

where the index a determines, as above, the spin compo- 
nent of the operators, a € {x,y,z}. The operators are 
represented in the interaction picture (see the Appendix), 
which means that their time-dependence is governed by 
Hq, i.e., S(r„,i) = e*^o*S(r„)e-*^o*. We assume that 



the nuclear spins evolve on a time scale which is much 
longer than that of the electron spin dynamics, so that 
for all times of interest we can use K„(i) = K„. We see 
that, if we focus on the dynamics of one specific nuclear 
spin (as we will do below) , our perturbative approach re- 
quires that Avo is much smaller than the energy scales 
in the unperturbed Hamiltonian Hq. 

We then substitute the perturbation ([5]) into the equa- 
tion of motion (|4]), make use of the fact that a trace is 
invariant under cyclic permutation of the operators in its 
argument, and then finally arrive at 



dp'' 



-{AvorY.{f(s''( 



(6) 



(r„,t)5^(r™,i'))cdt'K,</5l- - / (^''(r™,i')^'^(rn,i))cdi'K,/5'^<] 



where now brackets are used to denote the trace over the 
electronic part of the density matrix, (...) = Tre{. ■ ■ p^} 
in general, and (. . . )c = Tre{.../5g} for the fully con- 
nected terms. From ([6]) we can derive equations of motion 
for the expectation values of the nuclear spin operators. 

A single nucleus, K — 1/2 



with the vector S defined as S = Auo(S(r„, f)). Obvi- 
ously, the traces appearing in ([7|) and the corresponding 
averages (...) involve tracing and averaging over the nu- 
clear spin part of the system only. Since no confusion is 
possible, we do not label these traces and averages sepa- 
rately with a k. 



Let us first focus on the transparent case of one single 
nucleus with nuclear spin K = 1/2, interacting with a 
local electron spin density. Up to first order, i.e., using 
only the first term in (|6]), we find 

^"=^'{^'^4} = -^'^'<l*"-'^-'l-^') (7) 

describing a precession of the nuclear spin around the lo- 
cal electron spin density at the position r„ of the nucleus, 



The exchange of angular momentum between electrons 
and the nucleus is to leading order described by the sec- 
ond term in (jS]), i.e., second order in H^. We write 
the correlation functions of the electron spin in terms of 
the quantities R"''' and x"'' introduced in ([3]), and make 
again use of the cyclic invariance of the trace. For the 
convenient case of JiT = 1/2, we can simplify products of 
nuclear spin operators using K'^K^ = ^e°''"^K'^ + ^6°"^. 
Along these lines we simplify the second order term. 



dt 



(2) 



^^^"^^ {(fl''" + ix^'')TT{k''k''k''p'' - k''k''p''k''} - (i?'"= - ix'")Tr{i^"ii:V^' - k'^p'^k^k''}} 

[Avof 



k''] - , X"] - ) + zx'"( [[if" , i^l _ , k' 



(Avo)' 



^^R^'^ik'') - R^^k") + £"''=ix^'}, 



(8) 



which is, combined with Eq. ([7]), identical to the expression given in Eq. ([T]). 
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Many nuclei, K >l/2 

When following the same derivation for the case of many nuclei and spin higher than K — 1/2, it is generally not 
possible to derive a closed set of equations for d(K.n) /dt. The first order equations do not change and are still given 
by Eq. ([7]). The second order equation however does change, and becomes 

(2) 

d{K^ = ^-^{n'^ki) ~ R'\k'},) + e'^'^x'l{[Ktk^^^^^ (9) 



where the correlator x has acquired two more indices, 

Xt^^-^ j\[S''{Vn,t),S\v^X)]-)cdt'. (10) 

This correlation function now describes the non-local sus- 
ceptibility of the electron spin density, i.e., the linear re- 
sponse in the spin coordinate 5° at the position r„ due 
to a perturbation along at another position r,„. 

We notice two practical problems with Eq. First 
of all, for K > 1/2 we cannot use the relation k^K'^ = 
^^abcf^c _|_ igab simplify products of nuclear spin op- 
erators, and secondly, also correlations {k^K^) between 
different nuclei play a role. This makes it impossible to 
derive from Eq. dH) a closed set of equations for (i(K„) /dt. 

However, typically the number of nuclear spins in the 
ensemble is large {N ~ IC'-IO^). When focusing on 
the dynamics of one single nuclear spin, as we do in 
Eq. ([9]), the presence of the large number of other nu- 
clear spins can, to first approximation, be considered 
as a static 'classical' field. In other words, we can as- 
sume the fluctuations in the nuclear field {k^K^)c to 
be much smaller than the average values of the field 
{k^){K^). This assumption then allows us to replace 
(kik^) {kt){k^) in the last term of Eq. ©. We 
emphasize here that in a standard FGR calculation, the 
same assumption is implicitly used: Transition rates be- 
tween different spin states of a single nucleus are evalu- 
ated, ignoring any possible corrections due to correlations 
between this spin and the other nuclear spins. 

With this simplification, the last term of Eq. ^ 

now reads (^^^o)^ Em e"^''X^l(^m)(^«) and describes 
an effective electron mediated internuclear spin-spin cou- 
pling. This effect can be interpreted as a correction to 
the first-order equation ([T]): The effective field around 
which the n-th nuclear spin precesses is adjusted 5*" — 
S^ + iAvo)"^ Em Xnmikm), i-^., it acquires a contribution 
due to the effective field produced by all other nuclear 
spins, mediated by the electron spin 22 1. 

The third term in (|9]) contains products of spin op- 
erators of the same nucleus, and must thus be treated 
on different footing. In most realistic experimental cir- 
cumstances, the temperature is much larger than the 
nuclear Zeeman energy. This means that the (ther- 
mal) equilibrium density matrix of the nuclear spins 



can be approximated to be isotropic, i.e., /5q oc 1. As 
long as we focus on small deviations from this equilib- 
rium, i.e., small polarizations {K°'), the nuclear spin den- 
sity matrix is always close to isotropic. In this case, 
quadratic terms like {K^K^)c will dominate the contri- 
bution from (^[K^, K^]+)c in since they are of or- 
der unity, ((if°)^)c w lK{K + 1). We thus can con- 
centrate on these quadratic contributions, and approxi- 
mate (k^k^) Pi ^S'"'K{K+l). We see that then finally 
Eq. (0) reduces to the simple result of Eq. ([5]). 

Let us repeat here that this reduction is based on two 
assumptions: (i) we neglect the correlations between the 
nuclear spin under consideration and all other nuclear 
spins (as is also done in an FGR calculation) , and (ii) we 
only allow for small deviations of p'^ from equilibrium, 
so that {k^k^) « \5'"^K{K + 1). In fact, one does not 
necessarily have to use these assumptions and make this 
simplification. Equation ^ gives the time-evolution up 
to second order in ff^f of the full nuclear spin density ma- 
trix in terms of the electron spin fluctuations and suscep- 
tibilities. One could solve this Equation (numerically) for 
any initial condition /5q and calculate the time-dependent 
expectation value of any desired nuclear spin operator. 



IMPLEMENTATION 

With the formalism presented in Eqs ([I])-©, we are 
able to express the hyperfine driven nuclear spin dynam- 
ics in terms of the susceptibility and fluctuations of the 
electron spin density. Let us now explain how one could 
calculate the necessary electronic correlation functions 
using linear response theory. For a given setup, we write 
down the full set of Bloch equations describing the evolu- 
tion of the electronic density matrix under Hq , and solve 
for the steady state solution p^") = pg. The elements 
S"' are then simply found as the equilibrium expectation 
values 5" = Avo{S''{vn))pW ^ Tr{^'^(r„)p(o)}. Next, 
we add to the set of Bloch equations the first-order effect 
of a perturbation H' = A ■ S(r„), i.e., we add the terms 
i[/5'°\ and solve for the new steady state density 
matrix p'-^f The spin susceptibility x"'' needed in Eq. ([1]) 
is then found as the part of (S'°(r„))p(i) = Tr{S'°(r„)p(i)} 
which is linear in A**. The fluctuations i?"** are found in a 



similar way, the terms to add to the Bloch equations then 
read [/5(°\ 7?'] + -2Tr{i7'p(°)}p(°), where the last term re- 
moves the connected part of the electron spin correlators. 
We will illustrate the procedure below. 

From this short outline, it is already clear that this 
approach has several properties which are different from 
a traditional FGR calculation: (i) For an FGR calcula- 
tion, there should be an obvious basis to calculate all 
occupation probabilities in. Here we can choose any ba- 
sis which is most convenient for evaluating the steady 
state electronic density matrix, (ii) When calculating 
separate FGR transition rates, it should a priori be clear 
which quantization axis to choose for the nuclear spins. 
In our formalism, as soon as /S*-^-* is found for the two 
types of perturbations, all elements of x ^-nd R can be 
read off immediately, yielding the dynamics for the full 
vector (K„), not just the polarization along a single axis, 
(iii) For FGR to work, all levels involved in the hyperfine 
transitions (including possible virtual states in higher or- 
der processes) should be well separated in energy, i.e., the 
splittings should be larger than all decay rates and other 
possible incoherent processes present. In our method all 
incoherent processes are accounted for in the set of Bloch 
equations, we do not have to compare them with the 
other energy scales, (iv) In an FGR calculation, when 
focusing on the nuclear spin dynamics of a small subsys- 
tem being part of a larger system (e.g. a quantum dot 
coupled to multiple reservoirs), one usually includes only 
the energy splittings inside the subsystem in the actual 
perturbation theory. All other processes relevant for the 
electron spin dynamics (such as the coupling to the reser- 
voirs) have to be accounted for 'afterward' in an effective 
density of states, which is not always possible. In our 
approach, we treat the coherent and incoherent parts of 
the dynamics in the subsystem on equal footing. One ad- 
vantage is that the perturbation theory in our approach 
does not break down when the energy splittings in the 
subsystem vanish, as long as the energies corresponding 
to the incoherent dynamics are large enough. 

Let us mention here also a drawback of our method. 
Our main result ([l} describes only the dynamics of aver- 
age nuclear fields (K). Temporal fluctuations and other 
stochastic properties of the fields are not accessible via 
this approach. In an FGR calculation, one can separate 
the rates for flipping nuclear spins up and down, and 
derive from this separation a Fokker-Planck equation de- 
scribing the stochastics of the nuclear fields . 

EXAMPLE APPLICATIONS 

We will now illustrate the above with a calculation of 
the nuclear spin pumping in two example systems. First, 
we will consider nuclear spin pumping in a trivial toy 
model concerning a single quantum dot coupled to two 
leads. We will show how a smart implementation of FGR 
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can produce an identical result in this case. Then we 
will investigate a more complicated setup, involving a 
double quantum dot, which cannot be dealt with in an 
FGR approach. We will show however how our formalism 
can be applied without any problem, yielding a general 
equation for the nuclear spin pumping in this setup. 

Single dot: Fermi's golden rule with modifications 

We consider a single quantum dot connected to two 
leads, as illustrated in Fig.[TJ Due to the strong Coulomb 
repulsion, only one excess electron can occupy the dot. 
If we apply an external magnetic field Bext = Bg^tz, the 
energy levels for spin up and down in the dot are split by 
the electronic Zeeman energy. We take the tunnel rate 
into the dot Fjn to be very fast and equal for both spin 
directions. The outward tunnel rates F-i-^ are different 
for the two spin directions (to provide a physical picture 
with this assumption: one could assume the right lead 
to be ferromagnetic). Since the electron occupying the 
dot will thus have on average a non-zero spin polariza- 
tion along the z-direction, hyperfine interaction with the 
nuclear spins is expected to produce a non-zero nuclear 
spin polarization along this axis. 

In typical experiments, the temperature ('^ 100 mK ^ 
10 fieV) is much smaller than the orbital level spacing 
in the dot (^ 1 meV). In that case, the electron will 
in principle occupy the orbital ground state of the dot, 
which has an s-type character. As long as also the energy 
scale of the hyperfine interaction is much smaller than the 
orbital level spacing, we can project Hm to the orbital 
ground state, and write 

HM{t)^Y.^^^^"(*)^n, (11) 

n 

where the hyperfine coupling coefficients are defined as 
An = ^wo|i/'o(rn)P, i-C, in terms of the envelope wave 
function "00 (r) of the orbital ground state. The operator 
S is now the electron spin operator, instead of spin den- 
sity. A usual further simplification to make is to assume 
that |-0o(rn)P is approximately constant in the dot. This 
implies that |'(/'o(r«)P ~ (Nvq)-'^, where N - lO-^-lO^ is 
the number of nuclear spins in the dot. In this case, all 
coupling coefficients simply reduce to An = A/N . 

Let us now calculate the nuclear spin pumping rate 
along the z-direction, d{Kn) j dt, for this system. Since all 
nuclear spins are coupled to S with the same coefScient, 
(Kn) fo^' the n-th nucleus equals the ensemble averaged 
nuclear spin (K^)- We only take into account a finite 
polarization in the z-direction, so we use the equation 

(12) 

with the functions = AlK{K + l)(x''^ - X^"") and 
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FIG. 1. Energy diagram for a trivial single quantum dot setup 
in which nuclear spin pumping is to be expected. The dot is 
coupled to two leads with different chemical potentials. One 
single-electron level lies within the bias window, which is split 
into two sublevels by the Zeeman energy (we assumed a nega- 
tive p-factor). The tunnel rate from the left lead into the dot, 
Fin is large and equal for both spin directions. The outward 
tunnel rates F^^ are different for spin up and down. 



Q^^ = Al{R^'^ + Ryy). Note that x and R contain now 
correlators of the electron spin, not electron spin density. 

A finite nuclear polarization in the z-direction behaves 
on the time scale of the electronic dynamics as a static 
contribution to the magnetic field experienced by the 
electrons. We incorporate this contribution into an ef- 
fective field Bcff = i?oxt — A{K^-) (we assume a negative 
^-factor and a positive hyperfine coupling constant, as is 
the case in GaAs). The coherent time-evolution of the 
electronic 2x2 density matrix of this simple system is 
then given hy dtp — i [i?off 5*^ , /5] _ , where we express BcS 
in terms of a frequency. To this we add the incoherent 
tunnel rates F-i- and T^, which, to a good approximation, 
include the presence of the leads and the left and right 
tunnel barriers into the equations of motion for the elec- 
tron on the dot. Together, this yields the set of equations 



dtpti = {iBeS ~ ^'r)pn 
dtPn = {-iB^ff - ^T)pi^, 



(13) 



where F 



F^ + Fj^. These equations can be solved for the 



stationary situation dtp^^^ — 0, yielding p*-^'' — 

(0) 



and 



yO||.' — 0. We only need the response functions 

X'^y and x^"", so we add the terms A^S'^ + AySv]^ 

to Eq. (|13p and solve for the new stationary solution p^^\ 
The correlator x^^ is then simply given by the term in 
Tr{5^p(i)} linear in A^ , and in a similar way we find x^^ ■ 
The results are combined to 



and yielding 

R^^ 



Ryy 



2R2 I lp2 ■ 



(15) 



Combining all together using Eq. p^ . we thus find 



dt 



F^(i-(i^^))-F^(i + {if^)) 



4(Be 



(16) 



where for simplicity we assumed K = 1/2. 

Let us compare the result (|16p with the pumping rate 
we get using an FGR approach. We focus on transi- 
tions between the two states |t-U-n) and IJ^ftn), where the 
double arrow indicates the nuclear spin along the z-axis. 
Neglecting the nuclear Zeeman energy, the two states are 
separated by the energy fiBeS, and both are broadened 
by their respective electronic decay rate r^(^-). Hyper- 
fine induced transitions from It-IJ-n) to Wri^n) contribute 
positively to d{K^)/dt, and their rate is given by 

F+ = 2^1 iiU \^Aj-k+\ Un) ?P^V, (17) 

the square of the relevant matrix element, multiplied by 
the chance pi of finding the system initially in It-IJ-n) and 
the effective 'density of states' for this transition T). 
For the occupation probability pi, we take pi — [\ — 

(K^-)) X p'^K In this simple case, it is possible to incorpo- 
rate all incoherent dynamics (the two decay rates) into a 
smartly chosen V: We take the sum of the level broad- 
enings of the two electronic levels |t) and H). Assuming 
a Lorentzian shape for these broadenings, we thus use 



V 



so that we can write 



1 



iF 

2^ 



lp2 ' 
4^ 



f+ = a: 



2 ^iil 



F2 



(18) 



(19) 



We derive similarly an equation for F_ and then combine 
the two to write for the net nuclear spin pumping rate 



dt 



= A 



2 ^dh 



{k^))-T^{l + {k^)) 



A[B,^,-A{K^)Y+T^ 



(20) 



We see that this result coincides with Eq. (fT6| . The fact 
that the FGR approach in this case also works for small 
energy splittings, i.e., in the regime B^s ^ E, is due to 
the fact that all incoherent effects in this model can be 
incorporated relatively simply into V. 



X 



yx 



2S2^+iF2- 



(14) 



Now we evaluate R^^ and R^^. To this end, we add the 
terms [p^") , A"" 3"= + Av Sy]+ - 2Tr{ {A'= S"" + Ay Sy)p^°'> }p^°^ 
to Eq. (IT51) . and then look for the linear responses in S*^ 



Double dot: Beyond Fermi's golden rule 

Let us now illustrate what can happen when the setup 
becomes more complicated. We consider a double quan- 
tum dot connected to two leads with different chemical 
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FIG. 2. Energy diagram of a more complicated setup. Two 
quantum dots are coupled to each other and to two leads with 
different chemical potentials. The two dots have different ef- 
fective g-factors, so that the application of an external mag- 
netic field yields different Zeeman energies in the dots. The 
double dot is tuned to the (0,0) -> (f,0) (0, f) -5> (0,0) 
transport regime. The outward tunneling rates F are equal 
and much smaller than the inward rate Fin. 



potentials, as depicted in Fig. [51 The dots are tuned such 
that an electronic transport cycle involves the charge 
states (0,0) (1,0) (0,1) (0,0), where (n, m) de- 
notes a state with n{m) excess electrons on the left (right) 
dot. The electronic levels in the two dots are both split 
by a Zeeman energy in a spin up and spin down level. We 
assume that the splitting in the left dot is larger than in 
the right dot (e.g. due to a size-related difference in ef- 
fective g- factors), and that the (1,0) spin down level has 
a too high energy to play a role. The (1,0) and (0, 1) 
spin up levels are detuned by an energy d and coupled 
to each other with a tunnel coupling t. We again assume 
the outward tunneling rates to be much slower than the 



inward, F ^ Fi„, so that the system effectively will never 
occupy the (0, 0) charge state. 

We will focus on nuclear spin pumping in the right dot. 
A non-zero nuclear spin pumping rate is expected along 
the direction of the magnetic field since the average local 
electron spin polarization along this direction is finite. 
Denoting the level in the left dot by \L) and the two 
levels in the right dot by |t) and ||), the Hamiltonian for 
the three-level system reads 

Ha^d \L) {L\ + Beff li) (il + t \L) (;| + t |i) {L\ , (21) 

where the feedback of the nuclear spins on the electron 
spin dynamics again is incorporated into the effective 
magnetic field B^s = ^cxt — A{K^). Using this Hamilto- 
nian and adding the effect of the incoherent decay rate 
F, we can write down the equations of motion for p, and 
solve for p*^"' , yielding 



(0) 

Pl 



+ 4i2 + f2 



4d2 + 8t2 



(0) / 

Pit = 



. p2 ' 

Adt 



(0) 

p\ - 

2itr 



4^2 + 8t2 + r2 ' 



(22) 



4d2 + 8^2 + r2 ■ 



and zero for the other five elements. We then proceed 
as we did above, and add consecutively the perturba- 
tions z , iJ'] _ and [p(o),i?']+-2Tr{F'/5(o)}/5(°) to the 
equations of motion. We solve in both cases for the new 
steady-state solution p'^^-' and extract the linear responses 
needed. Combining all together, we find 



dt (4d2 + 8i2 + F2) [4(Beffd - B2ff + t2)2 + (552^ _ ^B,sd + 4d2 + 4i2)r2 + F*] ' ^ > 

where — B^xt ~ A(K^). This equation gives the nuclear spin pumping in the right dot along the z-axis, without 
any restriction imposed on the relative magnitude of the parameters B^ff, d, t, and F. 

Let us now try to evaluate the same pumping using an FGR approach. The most transparent electronic basis to use 
is {\L) , It) , li)} since in this basis we have well-defined incoherent in- and outward tunneling rates. We see however 
that the stationary density matrix in this basis (j22p has off-diagonal elements. In order to proceed, we have to make 
several assumptions. Let us consider the limit of F In this case the off-diagonal elements of can be neglected, 
and we can approximate w \L) {L\. We further assume that we have three well-separated levels, i.e., Bctf, d ^ T. 
In this limit the nuclear spin pumping is resulting from a second order transition: (i) tunneling from the initial state 
\L itn) to the virtual state |4,'|t„), and then (ii) an electron-nuclear spin fiip-fiop from \ii[n) to the final state It-D-n)- 
We thus use a standard second order Fermi golden rule to write for the pumping rate 

^ ^ iE,,.~E,,J^ (2^) 



where Ht = t \L) {W + t \],) {L\ is the tunnehng part of 
the Hamiltonian Hz- The occupation probability of the 
initial state is pL-i\„ = ^ — (-^^)- The 'density' T) for this 



transition is mainly set by the level broadening of the 
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final state, so V = ^T/[A{Bcs - df + T^]. This yields 



dt 



d2[4(Bcff-rf)2 + r2]' 



(25) 



which coincides with (|23p in the limit Bcff, d ^ T ^ t. 
(Note that Bcs — d F is not a necessary condition.) 

There are other limits in which considering the right 
rates to calculate, combined with a good approximation 
for p^^^ and a smart choice for D, can result in the right 
pumping rate. In general however, it is not trivial to un- 
derstand how to incorporate all incoherent dynamics into 
v. Even if one would simply diagonalize the stationary 
density matrix p^*^^ and try to evaluate all hyperfine tran- 
sitions in the new basis, one would have to transform the 
incoherent processes in correct effective densities of state. 
In the approach using the electron spin susceptibility and 
fluctuations, this is all done on the fly and its validity is 
not restricted to certain limits of the parameter space. 

Let us emphasize here that the fact that a correctly 
adapted FGR approach did work in our single dot ex- 
ample but not in the double dot example, is only due to 
the complexity of the coherent and incoherent electron 
spin dynamics in the two cases. The success of the FGR 
approach is not inherently connected to the number of 
quantum dots in a setup: Also for more complicated sin- 
gle dot systems, the FGR approach might break down. 

The two examples given above both consider localized 
electron spins. Our approach however, works in princi- 
ple equally well for systems with delocalized electrons. 
In this CStSG, cLS explained above, the susceptibility and 
fluctuations needed are correlators of electron spin den- 
sity instead of simply electron spin. This however does 
not change the complexity of the formalism: One only 
has to find the right set of (Bloch) equations to describe 
the electron spin density. 



CONCLUSIONS 

To conclude, we presented a new approach to evaluate 
the nuclear spin dynamics driven by hyperfine interaction 
with non-equilibrium electron spins. To describe the dy- 
namics up to second order in the hyperfine coupling, it 
suffices to evaluate the susceptibility and fiuctuations of 
the electron spin. This approach does not rely on a sep- 
aration of electronic energy scales or the specific choice 
of electronic basis states, thereby overcoming practical 
problems which may arise under certain circumstances 
when using a more traditional formalism based on rate 
equations. 
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APPENDIX: DERIVATION OF EQ. (gD 

Here we will show how to derive the time evolution 
equation ^ using second order perturbation theory in 
the interaction picture. The total Hamiltonian of the sys- 
tem H is split into an unperturbed part Hq and a small 
perturbation H'. In the interaction picture, all opera- 
tors acquire a time-dependence governed by Hq, so that 
for any operator A we have Ai{t) = g^^ot j^^-tHot 
assumed for simplicity a time-independent Hq). The re- 
maining part of the total Hamiltonian, the perturbation 
H', governs the time-dependence of the wave function, 
idt = Hj{t) IV'/), where Hj{t) is the Hamiltonian of 
the perturbation in the interaction picture. The time- 
evolution operator for \ipi) thus reads 



Ui{t,t')^T 



exp 



H'j{T)dT 



(26) 



where T is the time-ordering operator. The relation be- 
tween Uj and the 'regular' time-evolution operator in the 
Schrodinger picture U is straightforward to derive and 
reads J7(i,<') = e-'"°*Ui{t,t')e'"°*' . 

Let us now consider the evolution of the density ma- 
trix of a system described by the Hamiltonian Hq + H'. 
We assume that at some time in the past to ~> ~oo the 
system was still unperturbed and could be described by a 
stationary density matrix po which is determined solely 
by Hq. The perturbation is then switched on adiabat- 
ically, formally done by replacing H' e^^H' , where 
rj is an infinitesimally small positive real number. The 
density matrix at any later time t can then be written 
as p{t) = U{t,to)poU{to,t), where U describes the time- 
evolution due to Hq + e^^H' . In terms of evolution in the 
interaction picture, we see that we can write 



-iHnt 



°'Ui{t,to)pQUi{tQ,t')e 



iHot 



(27) 



where we used that e^^°'^°poe '-^"'f = po since po is per 
definition stationary under Hq. 

Introducing pi{t) = e^^"* p{t)e~'^^°* , the 'interaction' 
density matrix, we note that it obeys the simple relation 

Pi{t)^Ui{t,tQ)poUi{to,t'), (28) 

and that for the expectation value of any observable 

{A{t)) = Tr{Ap{t)} = Tr{Ai{t)pj{t)}. (29) 

From Eqs (I28|) and (|26|) it follows immediately that we 
can write as time-evolution equation for pi{t) 



dpi 
dt 



-i[H'j{t),pj{t)]^. 



(30) 
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If one knows the stationary density matrix po, one can 
iteratively use Eq. (j30[) to approximate 



dp 
~dt 



' =-^[H'J{t),po]- 



[H'J{t),[H'J{t'),po]-]-dt'■ 



m 



In our case, we would like to obtain a time-evolution 
equation for a reduced density matrix, i.e., separate the 
system into two subsystems, 5'-^-' and 5^^^, and trace out 
the degrees of freedom of one of the subsystems. We 
assume that initially, at to, the density matrix can be 
separated, po = pg^^ (8) p'i^\ where the superscripts label 



the subsystems. We then trace in Eq. (|3Tt over the de- 
grees of freedom of 5*^^-', which yields up to second order 
in H' for the reduced density matrix 



dt 



Tr2{[i7;(t),[i?;(t'),/5W®4'V]-}rft' 



In this expression, tracing out the degrees of freedom of 
5(2) general introduces correlation functions of those 
coordinates of 5*^^-' to which H' is coupled. If one sepa- 
rates the contributions of the connected and unconnected 
correlators thus appearing in this expansion, one can see 
that all not fully connected correlators are actually the 
higher order iteration terms for lower order perturbation 
terms, which can all be collected together by replacing 
Pq^ p^^^(t), similar to all terms in the iteration series 



([31]) combining to ((30|) . In terms of only fully connected 
correlators of the coordinates of 5^^^ we can thus write 



dp) 



(1) 



= -,:[Tr2{if;Wp^'^},P«W]- 



dt 



Tr- 



(33) 



indicates that only fully connected contributions should 
be included here. Note that we dropped all 'interaction' 
indices I in the main text: all operators and density ma- 
trices are implicitly written in the interaction picture. 
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which equals Eq. (U). The line over the trace symbol 



